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Abstract: Using a microscopic approach beyond RPA to treat anharmonici- 
ties, we mix two-phonon states among themselves and with one-phonon states. 
We also introduce non-linear terms in the external field. These non-linear terms 
and the anharmonicities are not taken into account in the "standard" multi- 
phonon picture. Within this framework we calculate Coulomb excitation of 
208 Pb and 40 Ca by a 208 Pb nucleus at 641 and lOOOMeV/A. We show with 
different examples the importance of the non-linearities and anharmonicities 
for the excitation cross section. We find an increase of 10% for 208 Pb and 20% 
for °Ca of the excitation cross section corresponding to the energy region of 
the double giant dipole resonance with respect to the "standard" calculation. 
We also find important effects in the low energy region. The predicted cross 
section in the DGDR region is found to be rather close to the experimental 
observation. 



1 Introduction 



States that can be interpreted as the first quanta of collective vibrations are a 
general property of quantum mesoscopic systems which can be found in various 
fields of physics. In nuclear physics, such vibrational states of the nucleus have 
been known for many years M. These one-phonon states are present both in 
the low-lying excitation spectra of nuclei and at higher energies. The latter are 
the Giant Resonances (GR). The existence of two-phonon states, i.e. states 
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which can be described as double excitations of elementary modes, has also 
been predicted since the early days of the collective model Q. Such states 
have been observed long time ago in the low-lying spectra. More recently, 
two-phonon states built with giant resonances have been populated in heavy 
ion inelastic scattering 0, in double charge exchange (n^,^) reactions 
and in Coulomb excitation at high energy For a review, see ref. [0. 

In the harmonic approximation, these states are predicted as degenerate mul- 
tiplets located at an excitation energy equal to the sum of the individual 
phonon energies. When the residual interaction is taken into account, the de- 
generacy is broken by the coupling between phonons. In the present article we 
consider the residual interaction of two-phonon states among themselves and 
with one-phonon states. Therefore, the eigenstates are linear combinations of 
one- and two-phonon components while the energies are shifted and splitted 
with respect to the harmonic limit. We will call such states mixed states. Ev- 
idence of such anharmonic behaviour can be found, for example, in a (7, 7/) 
experiment || where the observation of large dipole strength in the low-lying 
spectrum of some Sn isotopes is reported. Such strength is interpreted by the 
authors as due partially to the population of the 1~ member of the quintuplet 
of states based on the |2 + <g> 3~ > two-phonon state, and partially to the ad- 
mixture of the (one-phonon) GDR in the wavefunction of the state observed 
around 3.5 MeV excitation energy. As it has been shown in ||, the inclusion of 
the residual interaction among two-phonon states leads to small, but sizeable, 
anharmonicities also in the high-lying spectrum, namely for those states that 
in the harmonic limit are described as double excitations of GR. The mixing 
between one- and two-phonon states further increases the anharmonicities. 
When an external field acts on a nucleus, it excites the eigenstates of the 
internal hamiltonian, which, in our approach, are superpositions of one- and 
two-phonon states. 



The microscopic theory suited for the description of collective vibrational 
states is the Random Phase Approximation (RPA). Two-phonon states and 
their mixing among themselves and with one-phonon states can be generated 
by using boson mapping techniques and by taking into account terms of the 
residual interaction which do not enter at the RPA level In this way 

one has an RPA based approach to treat anharmonicities. 



In a nucleus-nucleus collision, the mutual excitation of the two partners is 
described as due to the action of the mean field of each nucleus on the other 
one, i.e. by a one body operator. Assuming that it induces small deformations 
of the density, only the particle-hole (ph) terms of the external mean field are 
usually taken into account. This amounts to consider as elementary processes 
only those corresponding to the creation or annihilation of one phonon. In 
this approximation, the external field is linear in the creation and annihilation 
operators of phonons. When the particle-particle (pp) and hole-hole (hh) terms 
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of the external field are also included, the direct excitation from the ground 
state to two-phonon states as well as the transition between one-phonon states 
become possible. These terms can be expressed as quadratic in the creation 
and annihilation operators of phonons and so correspond to non-linear terms 
in the excitation operator. 



In the "standard" approach, based on the independent multiphonon picture, 
the effects coming from both anharmonicities and non-linearities are neglected 
(see for instance ref.pi[|). Recent experimental data on Coulomb excitation at 
relativistic energies have raised some questions on the adequacy of that picture. 
Indeed, in the excitation of 136 Xe on 208 Pb, the experimental cross section to 
the double GDR (DGDR) has been found to be 2 to 4 times larger than the 
theoretical one Q. Recently, new experimental results || on the excitation 
of several nuclei have shown that the disagreement ranges from about 10% to 
60%, being about 30% in the case of 208 Pb. In a previous paper fll2] , by using 
a one-dimensional oscillator model to mimick nuclear states, we have shown 
that the effects of anharmonicities and non-linearities can lead to an important 
enhancement of the cross section in the energy range around twice that of the 
GDR. In this model neither spin nor parity were taken into account. Besides, 
only one type of phonons was considered. In the present paper we present more 
realistic calculations, where the collective states of the target nucleus and the 
action on it of the Coulomb field of the projectile are described starting from 
RPA. Anharmonicities and non-linearities are included by means of boson 
mapping techniques. Both low-lying collective states and giant resonances are 
considered as elementary phonons. 



We have done calculations for the 208 pb+ 208 Pb system at 641 and 1000 MeV 
per nucleon for which experimental data exist ||. We have also studied the 
Coulomb excitation of 40 Ca in the reaction 208 Pb+ 40 Ca at 1000 MeV/ A al- 
though there are no experimental data for this case. In both cases we consider 
as elementary modes all natural-parity RPA phonons whose multipolarity is 
lower than 4 and whose contribution to the associated energy weighted sum 
rule (EWSR) is larger than 5%. Then, we have built the residual interaction 
in the one and two-phonon space and we have diagonalized the hamiltonian 
in this subspace in order to define the mixed states \<f> a > . By solving the 
time dependent Schrodinger equation in this subspace we get the probability 
amplitudes for each of the \<p a > states from which we calculate the cross sec- 
tion. We will describe in detail the results for the 208 pb-|- 208 Pb system at 641 
MeV/A, the results at 1000 Mev/A being essentially the same except for the 
absolute values of the cross section which are higher in the latter case. We will 
mostly discuss the two regions around the energy of the states built with two 
low-lying states or with two giant resonances. We will see that non-linearities 
and anharmonicities may strongly change the cross sections associated with 
some specific states. Their influence is found to be of the order of 10% for 
208 Pb and 20% for 40 Ca in the DGDR energy region, bringing the theoretical 
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results closer to the experimental ones. 

In the next section we detail the model employed and in section 3 we describe 
the semiclassical electromagnetic field used to excite the nuclei. Section 4 is 
devoted to the description of the results on 208 Pb where we discuss in a detailed 
way the effects of both anharmonicities and non-linear terms on the excitation 
cross section. The results for 40 Ca are reported in section 5 and finally we draw 
our conclusions in section 6. 



2 The multiphonon picture 

Heavy ion collisions at high incident energies can be described within a semi- 
classical approach, where the relative motion is treated classically while quan- 
tum mechanics is used for the internal degrees of freedom of the colliding 
nuclei. For grazing and large impact parameter collisions the densities of the 
two nuclei have a small overlap. Therefore, the total hamiltonian can be writ- 
ten as 

H = H A + H B (1) 
where H A (H B ) denotes the hamiltonian of nucleus A(B) and 

H A = H° A + £ < a\U B (R(t))\a' > a) a a a , = H° A + W A {t) (2) 

aa' 

H A being the internal hamiltonian of A. W A describes the excitation of A 
by the mean field Ub of nucleus B, whose matrix elements depend on time 
through the relative coordinate R(t). The sums over the single particle states, 
denoted by a and a', run over both particle and hole states. 

2.1 Harmonic approximation 

Within RPA, the excited states > of each nucleus are described as super- 
positions of ph and hp configurations with respect to the ground state > 

\*v >=ql\^o >=E \ X ph^h ~ Y p u h a[a p ] | > (3) 

ph 

where the amplitudes X and Y are solutions of the RPA secular equation, with 
eigenvalues E v . The ground state is defined as the vacuum of the q u operators 

g,|^o>=o (4) 

In order to avoid unnecessarily complicated expressions, we do not introduce 
explicitly the coupling to total angular momentum and isotopic spin. When 
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the RPA phonons are mapped onto bosons [|13[], the internal hamiltonian of 
the nucleus can be written as 



H° = £ E V Q],Q V (5) 

V 

which shows that the excitation spectrum is harmonic. The boson operators 
Ql and Q v in the above equation are given by 

Ql = Y,i X ph B lh - Y ph B Ph] (6) 

ph 

with the same X and Y amplitudes as in eq.(3) but in terms of the boson 
images (B\ h and B p h) of the ph operators 

a\a h -> B ] ph + (7) 

In the above equation we have indicated only the first term of the boson 
mapping. Assuming that the external field induces only small deformations of 
the density, only ph and hp terms contribute to Wa and one gets 

W A (t) = E < P\U B (&(t))\h > a\a h + h.c. (8) 

ph 

By introducing the boson mapping of eq.(7), it can be rewritten as 

W A (t)=J2wl°(t)Ql + h.c. (9) 

V 

with 

=<y„\w(t)\v > (10) 

The independent multiphonon picture is based on eqs. (5) and (9). Within 
this picture, the Schrodinger equation can be solved exactly and the state of 
each nucleus at time t is found to be the coherent state 

|$(t) >= TT e -HM*)l 2 y ^j n " e- in " E ^(Qt) n "\0 > (11) 

v n v ' L v 

with 

t 

W W ^ e -in„E„t> dt , (12) 

where the integral is performed along the relative motion trajectory corre- 
sponding to a definite impact parameter. The probability amplitude to excite 
one- or two-phonon states is calculated by projecting eq. (11) on the corre- 
sponding states 

\u>=Ql\0> (13) 

K >=(i + <wr 1/2 QtQt'|o> (14) 
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where |0 > is the vacuum of the Q\, operators. Finally, the cross sections are 
obtained by integrating the relevant probabilities over the impact parameter. 

2.2 Non-linear excitation 

In this section we present an approach to go beyond the independent multi- 
phonon picture by eliminating its two main limitations, the first in the external 
field and the second in the internal hamiltonian. Let us first consider the pp 
and hh contributions to the sums in eq.(2), which are neglected in the lin- 
ear approximation for the external field. Assuming the same boson mapping, 
truncated at the lowest order, it is easily shown || that the mappings 

<{v-»E^ ( 15 ) 

h 

a h,a{' ^J2 B U B Ph' ( 16 ) 

V 

are exact, in the sense that they preserve the commutation relations between 
fermion-pair operators. Using these relations, the inclusion of the pp and hh 
terms in eq.(2) gives a W quadratic in the boson operators B ph and B ph . By 
expressing the latter in terms of the collective bosons Q} v and Q v one gets 

W = W m + E Wl°Ql + h.c. + E Wll,QlQ v , + E WZQlQl + h.c. (17) 

" vv' vv' 

where 

wl° = Y.(w vh x;* h + w hp v£) (is) 

ph 

is the standard linear response expression, whereas 

php'h' 

W%= E (Wrf8hh>-W h , h 5rt)x£Yfi, (20) 

php'h' 

provide new excitation routes. 

The matrix elements of Ub depend on the considered excitation mechanism. 
Since the general discussion we present here is independent of their form, we 
postpone to the next section the derivation of their expressions in the case of 
Coulomb excitation at relativistic energy. 

The hamiltonian Ha, with the inclusion of the terms W 11 and W 20 , is a 
quadratic form in the Ql and Q u operators. Therefore, a coherent state solu- 
tion to the Schrodinger equation still exists. We do not exploit this property 
because it does not hold any more when the anharmonicities are included, as 
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we are going to do in next subsection. The effects of introducing non-linear 
terms in the external field can be important, for example, whenever some se- 
lection rule disfavours one of the two steps necessary to make the transition 
from the ground state to a two-phonon state through the action of W 10 alone. 
The term W„„t describes the transition from the one-phonon state > to 
\ty u i > or from a two-phonon state to another one. In ref. fll4j-|l6l1 it was shown 
that these non-linear terms can lead to an increase of the population of two- 
phonon states. The term W 20 induces a direct transition from the ground 
state to a two-phonon state that can be very important. This effect has been 
already reported in || where the direct matrix element between the ground 
state and the dipole member of the low-lying |2 + ® 3~ > quintuplet of states 
in some Sn isotopes was found to be very large. Similar results, but involving 



double GR states, have been obtained in [17 



2.3 Anharmonic spectrum 



Let us now turn our attention to the other limitation of the independent mul- 
tiphonon picture we have stressed above, namely the assumption that the 
internal hamiltonian has the harmonic form of eq.(5). The simplest way to go 
beyond this approximation starts from the observation that in RPA only the 
phph and pphh terms of the residual interaction are taken into account. The 
pppp and hhhh terms, when expressed by the same boson mapping used be- 
fore, introduce a coupling between two-phonon states while the remaining, 
ppph and hhhp, terms mix one- and two-phonon states. Finally, when con- 
sidering two-phonon states one should also take care of the possible violation 
of the Pauli principle. In the boson mapping method the exclusion principle 
is introduced through high order terms in the boson expansion || built to 
conserve the fermion-pair commutation algebra. In such a way an additional 
residual interaction between two-phonon states coming from the particle-hole 
matrix elements is generated [plj. As a result of these different couplings, the 



eigenstates of the internal hamiltonian of each nucleus are 

l*« >= E c > > + E <, 2 h^ > (21) 

V U\V2 

Therefore, the states excited by the external field will be such mixed states and 
one cannot speak of pure one- or two-phonon excitations any more. However, 
when a |$ a > state has a strong overlap with a one-phonon state we may 
discuss the associated cross-section as part of the one-phonon cross-section. 
Conversely, if the |$ Q > state is dominated by its two-phonon components we 
may speak about two-phonon strength. 

For example, let us consider a state \<& a > which strongly overlaps with a 
two-phonon state, i.e. whose largest component is \uiu 2 >■ In addition to the 
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possible excitation of |$ a > via this two-phonon component, this state can also 
be excited by W 10 through its one-phonon components. However, the energy 
E a of |$ a > will be not far from E Ul + E U2 . Therefore, it will contribute to 
the cross section at that energy. In this sense, because of its structure and of 
its energy, one may say that it contributes to the two-phonon cross section. 



This fact was disregarded in ref. ||18|| , where the mixing of a huge number of 
one- and two- phonon states was considered. A good description of the width 
of the GDR was thus achieved. However, all the |$ a > states were considered 
to be one-phonon states when they were excited through their one-phonon 
components while the two-phonon excitations were calculated in |18|] as the 
transition to states of the form |$ Q (g) $ a / >. Therefore, that calculation is 
somewhat equivalent to consider a harmonic spectrum with the states |$q, > 
as the elementary quanta. 



2.4 Time- dependent excitation process 



The cross section is calculated, non perturbatively, by solving the Schrodinger 
equation in the space of the ground state and the I^q, > states. Then the time 
dependent state, \^(t) >, of the nucleus can be expressed as 

|*(0>=£A*(*) e ~ iE i$"> (22) 

a 

where the ground state is also included in the sum as the term a = 0. The 
amplitudes A a (t) are solutions of the set of linear differential equations 

Mt) = -iJ2e i{Ea - Ea,)t < ^a\W(t)\^ > (23) 

aa' 

and the probability of exciting the internal state |$ a > is given by 

P a = K(t = +oo)| 2 (24) 

for each impact parameter. Finally, by integrating P a over the impact param- 
eter we obtain the cross section 

+00 

a a = 2tt J P a (b)T(b)bdb, (25) 


where the transmission coefficient T(b) has been taken equal to a sharp cut- 
off function 6{b — b m i n ). The parameter b min is usually chosen such that the 
contribution from the nuclear part can be neglected. 
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3 Relativistic Coulomb excitation 



Let us now look in detail the multipole expansion of the external field of 
eq.(2). Alder and Winther have worked out an analytic expression for 



the Fourier transform of the semiclassical electromagnetic field in relativistic 
nucleus-nucleus collisions, with the assumptions that the projectile follows 
a straight-line trajectory and that the charge densities of both nuclei do not 
overlap. Therefore, to get the time dependence of the electromagnetic coupling 
potential the inverse Fourier transform of the expressions derived in |l9j can 



be taken. This procedure has the advantage that the multipole expansion of 
the time dependent coupling potential is readily known as well as its electric 
and magnetic components. 

Let us introduce the Fourier components of the time dependent coupling po- 
tential 

+00 

W(t) = — I e~ iut W(u)duj (26) 

— 00 

Introducing the expansion of the external field in multipoles W Xl1 

+00 

W(t) = —Y i [ (e iujt (-l) x+ » + e- tujt )W x ^(uj)du; (27) 
2n > J 

where we have taken into account the behaviour of the multipoles W Xfl for 
negative cu. 



It is shown in ref. JL9| that the contribution to W(u) of the (A,/i) multipole 
can be expressed in terms of electric (tt = E) and magnetic (tt = M) one-body 
operators 

W»(\u>\) = ^E^A M (-)(-ir^(M)v / 2ATT(-) A ^*(7r^A- Ai ) (28) 
t>7 ^7 v c 

where (3uj is the adiabaticity parameter related to the impact parameter b and 
to the Lorentz contraction factor, 7, and where are the modified Bessel 



functions. The expressions of the functions G^x^ can be found in ref. [19 
The 2 A -pole electric transition operator is given by 

M ( EuX ti = l^A+f) / JW ' V A L (^(^)^(r))rf 3 r (29) 
where J is the current density operator while j\ is a spherical Bessel function. 
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This operator can also be written down |13 

(2A + l)!!c A 



as 



M ( EooX/j) 



oo x (X + l) 



x ||p( r )r v (f)|-(rj A (^))+^J(r).rF A ,(OjA(^)}rf 3 r (30) 

where the charge density operator p(r) has been introduced. The second term 
will be neglected since, relative to the first, it is of the order of huj/2m p c 2 . 



To get the time dependent coupling potential as the Fourier transform (27- 
28) , we need the transition operators at any value of to. Since the argument 
of the spherical Bessel function in the operator is ur/c, the dependence on u 
and r of the multipole W Xfl will not factorize. Therefore there would be no 
factorization of the time and r dependence in the coupling potential. However, 
in the limit of long wavelengths, the first term in expression (|30D reduces to 
the well-known static electric multipole operator, 

M(EuXfi) ~Qx, = j p(r)r x Y Xfl (f)dr (31) 

which does not depend on u. In a similar way, the general magnetic operator 

M(Muj\fi) = - Z 2 ^ 1 ^^ 1 / J(r) • L(j x (^)Y x ,(r))dv (32) 

in the limit of long wavelengths, becomes 

M(MuXfi) ~ M x , = 1 J(r A J(r)) ■ (Vr x Y Xfl (f))dr (33) 

That means that, in the limit of long wavelengths, neither the electric nor the 
magnetic transition operators depend on u (we will therefore omit to in the 
arguments of M) and they will come out of the integral in (|27|) . We just need 
to know 

+oo 

Hxp(J3,t) = J (e^(-l) A+ ^ + e-^)u x K^u)du (34) 
o 

with n > 0. H XfM (j3,t) is an analytic function whose explicit expression can be 
found in the appendix. 

Therefore, in the long wavelength limit, the following analytic expression of 
the time dependent coupling potential 

mt) = ^£ G ^(-)(-ir^#A,(/M) M^x-ti) 

Z71V 7 V C A 

= E^(At)(-ir-M(7rA-/i)/e (35) 
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can be explicitly derived. From the use of the static multipole operators has 
followed that any of the term in the sum factorize into two elements, the first 
one depending on collisions properties, the second one acting on the nucleus 
being excited. 

If one considers only the electric components, which are the most important, 
only natural parity states can be excited in a first order calculation. However, 
in a coupled channel calculation, as the present one, non natural parity states 
have to be included since they can be reached, for example, through a two step 
process. However we will see in the following that these contributions remain 
small. 

We want to have a feeling about the limits that the use of the static electric 
operator imposes on the interpretation of our results. If a first-order harmonic 
and linear calculation was to be done, we could just compare the complete ma- 
trix elements < If\\A4(EuifX)\\Ii > and the approximate one < Jj||(5a||-^ >• 
The matrix elements connecting the ground state with one phonon states are 
consistent within a maximum of a few per thousand for low-lying states, and 
differ by a few percent when the GDR or the ISGQR are considered. 




O.O 20. 40. 60.0 SO.O lOO.O 

ct (fm> 



Fig. 1. Matrix element of the coupling potential between the ground state of 08 Pb 
and its GDR with magnetic quantum number zero, as function of time. There are 
two groups of lines corresponding to two different impact parameters. The solid 
lines have been obtained using the general expression (see eq. 30) of the electric 
dipole operator, while to get the dashed line the static expression (eq. 31) has been 
used. 

In a coupled-channel calculation not just a fixed Uif, but the full range of uj 
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40.0 



\ 




O.O 20. -40. SO.O SO.O lOO.O 

ct (fm) 



Fig. 2. As figure |T|, but for magnetic quantum number one. 



values will contribute to the Fourier transform ( p6| - p8[ ). As an illustration, let 
us consider the colliding system 208 Pb+ 208 Pb at Ei a b= 641 MeV to compare 
the exact multipole expansion with the long wavelength limit. The associated 
time-dependent transition matrices from the ground state to the giant dipole 
resonance, < GDR, fi\ W(t)\ >, are presented in figure [I] for the magnetic 
quantum number fi = and in figure |2| for \i — \ for two different values 
of the impact parameter. The solid lines correspond to calculations in which 
the general expression of the electric multipole operator has been taken into 
account, and the inverse Fourier transform of the corresponding amplitude 
has been carried on numerically. The dashed lines correspond to the use of 
the static electric operator and the analytic expression (|35|). We can see that 
qualitatively the time dependence is well reproduced, while the quantitative 
agreement gets better as the impact parameters increases. That is essentially 
due to the adiabatic cutoff that the modified Bessel function K^fiio) intro- 
duces. This function decays exponentially when the argument becomes bigger 
than 2 fl2TJ. Therefore, if the impact parameter increases the relevant range of 



u in the integral is reduced and we get closer to the long wavelength limit ex- 
pression for the interacting potential. This effect can be seen in figure ||| where 
the behaviour with the impact parameter of < GDR,ji = l\W(t = 0) 1 > is 
presented at t=0. This is the time at which the difference between both ap- 
proaches is maximum when A + \x is even. A similar behaviour is found when 
the matrix elements W 11 or W 20 are considered. 

The conclusion of this study is that, in the excitation energy region we will 
consider, it is reasonable to use the static electric operator. This amounts to a 
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35. 
b (fm) 



Fig. 3. Impact parameter dependence of the matrix element of the coupling potential 
at time t=0 between the ground state of 208 Pb and its GDR with magnetic quantum 
number one. 

considerable saving of calculations since, for each multipole, the time and the 
r dependence factorize. 



4 Results about the excitation of 208 Pb 



Let us now apply the above formalism to a specific nucleus, namely the Pb 
excited in a collision with a Pb nucleus at 641 MeV per nucleon. We will first 
discuss the effect of the anharmonicities on the RPA spectrum. Secondly, we 
will look at the effect of non-linear terms in the external field. Finally, we will 
consider the influence of both these terms on the excitation probability and 
the cross sections. 



4-1 Energy spectrum 



The one-phonon basis is calculated in the self-consistent RPA with SGII 
Skyrme interaction | 21| . Although we are using an explicit neutron proton 



representation the isospin results to be a rather good quantum number as far 
as collective states are concerned. We have selected all the states which exhaust 
at least 5% of the appropriate EWSR and, for a particular spin and parity 
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Table 1 

One-phonon basis for the nucleus 208 Pb. For each state its spin and parity, isospin, 
energy and percentage of the EWSR are reported. 



Phonons 



GMRi 
GMR 2 



GDRi 
GDR 2 



2+ 

ISGQR 
IVGQR 



3~ 

HEOR 



J7T J, 



E(MeV) 



\EWSR 



0+ 
0+ 



13.610 
15.022 



12.435 
16.662 



2+ 
2+ 
2+ 1 



5.545 
11.599 
21.815 



3.464 
21.302 



61 

28 



63 
17 



15 
76 

45 



21 
37 



(and isospin), we have grouped together the ones which are closer in energy 
according to the method described in ref. ||. We have considered the one- 
phonon states reported in table [I], i.e. the various components of the isoscalar 
monopole resonance (GMR), the components of the isovector dipole resonance 
(GDR), the low-lying 2 + state and the quadrupole resonances, both isoscalar 
(ISGQR) and isovector (IVGQR), and finally the collective low- lying (3~) and 
high-lying (HEOR) isoscalar octupole states. 

We have then constructed the residual interaction between the one- and two- 
phonon states and also among the two-phonon states. The two-phonon states 
are coupled to a total angular momentum and parity. In the case of the 1~ 
states, while the coupling between one- and two-phonon states is of the order of 
1/2 MeV up to 1 MeV, the coupling between two-phonon states is, in average, 
about one order of magnitude smaller. 

Then for each spin and parity the total matrix has been diagonalised in or- 
der to get the states |$ a >. Since these states are always dominated by one 
component we have decided to label them by the name of this dominant com- 
ponent. Table |2] gives for all 1~ states the total shift AE (in KeV) from the 
unperturbed energy Eq and their component on the GDR. Tables [3] and f| con- 
tain some information on the results of the diagonalization for the low lying 
and the high lying two-phonon states, respectively (see caption). We have re- 
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Table 2 

Characteristics of the \4> a > dipole 1~ states resulting from the diagonalization 
of the internal hamiltonian. In the first column we indicate the dominant com- 
ponent. The values in the second column remind us the energies associated with 
this component in the harmonic approach. The shift in the energy produced by 
the anharmonicities is indicated by AE (in KeV). We can compare these values 
with the diagonal matrix elements of the residual interaction, AEq (in KeV). In the 
last columns we report the amplitude with which the GDR's appear in such mixed 
states. 



Dipole 



States 



£ (MeV) AE (AEq) 



'GDK 2 



GDRi 






12.435 


-132. 


CO ) 


0.993 


-0.006 


GDR 2 






16.662 


—56. 


CO ) 


0.002 


0.994 


3~ 


<8> 


2+ 


9.009 


195. 


(200.) 


0.023 


0.000 


3- 




ISGQR 


15.062 


75. 


(67.) 


0.045 


0.000 


GDRi 


® 


2+ 


17.981 


-207. 


(-220.) 


0.043 


0.082 


GDR 2 




2+ 


22.207 


-23. 


(-36.) 


0.007 


0.048 


GDRi 


<8> 


ISGQR 


24.034 


33. 


(-10.) 


0.057 


-0.004 


3- 


<8> 


IVGQR 


25.278 


6. 


(4.) 


-0.014 


0.000 


GDRi 


<8> 


GMRi 


26.046 


18. 


(-27.) 


0.057 


0.000 


HEOR 


<8> 


2+ 


26.847 


25. 


(24.) 


-0.004 


0.000 


GDRi 


<8> 


GMR 2 


27.458 


-10. 


(-35.) 


0.039 


0.000 


GDR 2 


<8> 


ISGQR 


28.261 


-88. 


(-51.) 


0.007 


0.054 


HEOR 


<8> 


ISGQR 


32.901 


-31. 


(-30.) 


-0.004 


0.000 


GDRi 


<8> 


IVGQR 


34.250 


-44. 


(-47.) 


-0.007 


0.000 


GDR 2 




IVGQR 


38.477 


-174. 


(-174.) 


0.006 


0.000 


HEOR 


® 


IVGQR 


43.117 


-49. 


(-53.) 


-0.011 


0.000 



stricted these tables to natural parity states since the non natural parity states 
are essentially not mixed and weakly excited. Moreover, states with angular 
momentum greater than 3 have not been included since they do not play an 
important role in Coulomb excitation processes. From these tables one can see 
that the anharmonicities predicted by our microscopic calculations are small, 
the typical shifts in energy (AE) being a few hundred keV. Each multiplet 
appears to be splitted with a characteristic spreading equal to the global shift. 
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Table 3 

As table ^, but for the low lying states with natural parity. In the first column 
we give the dominant component while in the last one we report the second most 
important component and its coefficient. 



States 



E (MeV) 



AE (AE ) c conf Config. 



3.464 3" -256. (0.) 0.093 



5.545 



-364. (0.) 



-0.201 



6.927 



0+ 
2+ 



958. 
381. 



;il37.) 
(400.) 



-0.163 GMR 1 
0.195 2+ 



3- ® 2 H 



9.009 



195. 
161. 



(200.) 
(112.) 



-0.024 GDR X 
-0.091 3" 



2+ 



11.090 



136. 
178. 



(145.) 
(30.) 



-0.055 
-0.158 



3" 
2^ 



GDR ± 



17.981 



-207. (-220.) 
-4. (-4.) 



-0.083 GDR 2 
-0.010 GMRi 



The mixing coefficients are in average also small, around 0.05 and at maximum 
around 0.2. 



4-2 Excitation Processes 



Let us now study the characteristics of the excitation strength. We have seen 
that the excitation operator contains 3 parts. The first one is the linear re- 
sponse which is usually taken into account in the standard calculations: i.e. 
in the harmonic and linear picture. The strength associated to the operator 
W 10 is, in this picture, concentrated in the one-phonon states. The introduc- 
tion of a mixing between states with different numbers of phonons spreads the 
strength over more states. For instance, in the case of the GDR the strength 
will be distributed among the dipole states of tables [| in a fashion proportional 
to the c's coefficients. Analogously, the strength W 20 initially located around 
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Table 4 

Same as table |3[ but for mixed states with natural parity and with energies between 
22 and 29 MeV. 



States 



E (MeV) 



.r 



AE (AE 



^conf 



Config. 



GDR 2 


8> 2+ 


22.207 


1 


-23. 


(-36.) 


-0.046 


GDR 2 












3- 


—66. 




—0.018 


GDR 2 


Vy 




ISGQR 


(8) ISGQR 


23.198 


+ 


4. 


(3.) 


0.014 


GDRi 


<s> 


GDRi 








2+ 


35. 


(-15 ) 


—0.061 


ISGGR 






GDRi 


8) ISGG R 


24.034 


1~ 


33. 


(-10 ) 


—0.057 


GDRi 












3 _ 


-2. 


(-2.) 


0.010 


3 _ 


8> 


IVGQR 


3- 


8) HEOR 


24.766 


0+ 


34. 


(14 ) 


0.133 


GDRi 


8) 


GDRi 








2+ 


22. 


(-2.) 


0.225 


GDRi 


8> 


GDRi 


GDRi 


GDR\ 


24.871 


0+ 


41. 


(33.) 


-0.132 


3- 


8> 


HEOR 








2+ 


-189. 


(-192.) 


-0.223 


3 _ 


8) 


HEOR 


GMRi 


8> ISGQR 


25.210 


2+ 


42. 


(11.) 


-0.048 


ISGQR 






IVGQR 


8> 3- 


25.278 


l - 


6. 


(4.) 


0.015 


GDRi 












3- 


-24. 


(-25.) 


-0.018 


HEOR 






GMRi 


8> GDRi 


26.046 


1- 


18. 


(-27.) 


-0.057 


GDRi 






GMR 2 


8> ISGQR 


26.621 


2+ 


25. (8.) 


-0.033 


ISGQR 






2+ 


8> #£0/2 


26.847 


1" 


25. 


(24.) 


0.007 


GMR 2 


8) 


GDRi 








3 _ 


-39. 


(-44.) 


-0.018 


HEOR 






GMRi 


8) GMRi 


27.221 


0+ 


297. 


(56.) 


-0.128 


GMRi 






2+ 


8) IVGQR 


27.360 


0+ 


-196. 


(-195.) 


-0.016 


ISGQR 


8) 


IVGQR 








2+ 


-84. 


(-90.) 


-0.032 


IVGQR 






GMRi 


8 GDRi 


27.458 


l - 


-10. 


(-35.) 


-0.042 


GDRi 






GDR 2 


8> ISGQR 


28.261 


l - 


88. 


(51.) 


-0.055 


GDR 2 












317 


-60. 


(-62.) 


0.010 


2+ 





GDR 2 



GMRi 8> GMR 2 



28.633 



254. 



(74.) 



-0.100 GMR 2 8> GMR 2 



the two-phonon states, after the diagonalisation will be distributed over many 
states. 

Moreover, the various states have now two paths to be excited in one step, 
either through the W 10 excitation of their one-phonon component or via the 
W 20 interaction exciting directly their two-phonon part. Now, depending on 
the respective sign of the mixing coefficients, these two contributions may 
interfere constructively or destructively. 

In addition to these direct transitions from the ground-state, the term W u of 
the external field may induce transitions between excited states. These new 
excitation routes may modify the distribution of the excitation probabilities 
associated with different states. In the next subsection we will give a few 
examples where we will show the importance of the W 11 and W 20 terms and 
of the anharmonicities. 

4-3 Excitation cross-sections 

Let us now put all these ingredients together in order to compute the excitation 
probabilities and cross sections. All the natural parity states with angular 
momentum less or equal to 3 have been included in the calculations while for 
the non natural parity states we have included only the 1 + and the 2~ ones. 
By solving the coupled equations ([23|) we get the probability amplitude for 
each \<p a > state, from which we calculate the cross section by integrating over 
the various impact parameters associated with Coulomb inelastic excitations. 
The b min has been chosen according to the systematics of ref |2^| . We will 
describe in detail the results for the 208 Pb+ 208 Pb system at 641 MeV/A and 
we will first focus our discussion on the excitation of dipole states. 

In fig. [| we present the dipole excitation cross-section as predicted using var- 
ious approximations in order to disentangle the effects of the anharmonicities 
and non-linearities coming from W 11 and W 20 . We have run several calcula- 
tions corresponding to the various cases we can have, by switching on and 
off the different terms of the external field. From the figure it is clear that 
the spectrum is dominated by the dipole resonance. However, one can observe 
important modifications of the dipole strength for the different calculations 
compared with the harmonic and linear prediction. 

In particular, states which were not excited in the harmonic and linear picture 
can reach a sizeable cross-section when all the different corrections are taken 
into account. For instance, this is the case for the state around 9 MeV, which is 
mostly built out of the 1" component of the states resulting from the coupling 
of the low-lying 3~ and 2 + . In the first line of table |5|, the Coulomb inelastic 
cross-sections for this state at several degrees of approximation are given. One 
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208 Pb + 208 Pb, E lab = 641 MeV/A, L = 1 



Harmonic 



, 100 
10 

1 

100 
10 

! i 

b 

100 
10 
1 

, 100 
10 



b) 



c) 



d) 



W n =0, W 20 =0 I e) 



W"#0,W=0 EE f) 



w"=o,w 2 Vo i g) 



w"*o,w 2 Vo 



Anharmonic 



w"=o,w 20 =o 



w"*o,w 20 =o ! 



w"=o,w 2 Vo 



w"*o, w 2 Vo 



10 20 
E* (MeV) 







10 20 30 

E* (MeV) 



Fig. 4. Relativistic Coulomb target dipole excitation cross section for the 
208 Pb+ 208 Pb system at 641 MeV/A. Each bar corresponds to the cross section 
of a single state. 
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|2* <a> 3*= 



W (3) 
2* 



W (2) 



2 * w' 




W 1 °(2) 




g.s. 



1 2 tS> 3 > 

W 1<J (2) 

3" 



|2 <S> 3 



W*°(1) 



b) 



c) 



Fig. 5. Schematic representation of the Coulomb excitation of the |2 + ® 3 > state. 

can see that this two-phonon state is almost not excited in the harmonic and 
linear picture. Indeed, at this level of approximation, the most direct way to 
excite this state requires one E3 and one E2 transitions (see figure ^.a) which 
are not favourable. In this case the W n term does not help much because either 
we reach the state by one El plus two E2 transitions, as in figure |5].b, or by one 
E3 plus two El if in the first step we excite the 3~ state. In any case, at least 
one of the involved transitions is of high multipolarity. Conversely, the direct 
transitions due to the W 20 terms (see fig. |5].c) increases the cross section by a 
huge factor, bigger than 500. Indeed, this term is now a dipole transition which 
is strongly favoured. The importance of W 20 will decrease as the excitation 
energy of the state increases. For instance, the enhancement factor 500 reduces 
to about 50 for the dipole states \2 + ®HEOR > or \ISGQR®HEOR > whose 
energies are around 30 MeV. 

When the mixing of one- and two-phonons states is taken into account this 
state can be also populated by W 10 through its small GDR component (see 
table 0). In fact, although the c coefficient of the GDR component is small, this 
component gives a considerable contribution due to the fact that it is a one step 
dipole excitation. Moreover, the energy of the state (about 9 MeV) is lower 
than the one of the GDR state. All together the effect of the anharmonicities 
on the inelastic cross section is a factor about 100 times bigger with respect 
to our reference calculation. 

Finally, when all these different contributions are taken into account this dipole 
two-phonon state built from low-lying 3~ and 2 + is receiving 30 mb cross 
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Table 5 

Coulomb inelastic target excitation cross sections (in mb) for the 208 pp+ 208 Pb sys- 
tem at 641 MeV/A and for the mixed states which are identified by their dominant 
component (first column) and their angular momenta and parity (second column). 
In the third column is shown the reference result corresponding to a harmonic and 
linear calculation. In the fourth column the additional inclusion of only the W n 
non-linear term is allowed. Similarly, in the fifth column the only difference with 
the reference calculation is due to the addition of only the W 20 non-linear term. In 
the sixth column the results of an anharmonic and linear calculation are presented. 
The last column correspond to results of the anharmonic and non-linear approach. 



States 



harm. 
& lin. 



W 



11 



W 



20 



anharm. 



anharm. 
& non-lin. 



0.03 0.04 16.21 



2.60 



29.53 



ISGQR 



0.05 0.07 17.22 



3.63 



5.18 



22 < E < 28 (MeV) 1" 3.55 5.95 5.07 6.42 



12.18 



GDRi 



1.24 



2.07 



0.99 



7.64 



9.83 



ISGQR 



298.91 



332.56 



300.09 



278.35 



314.18 



section, while in the harmonic and linear limit it was just 0.03 mb. 

In this case the effects of non-linearities and anharmonicities interfere con- 
structively. That is not a general property. An example in which these effects 
interfere destructively is shown in table [5], where the excitation cross section 
to the dipole state \ISGQR<g>3~ > is given. In order to clarify this mechanism 
we have done a parametric calculation in which only three single phonon states 
were considered, namely the 3~, the 2 + and the GDR. Then we have mixed 
the single phonon \GDR > with the two phonon state |2 + <g> 3~ >, coupled to 
a total spin 1, in the following way 



|$i > = cos/?|2 + <g> 3" > + sin (3\GDR > 
|$ 2 > = -sin/3|2 + ®3- > + cos (3\GDR > 



(36) 



Increasing the parameter \f3\ we can go from a pure harmonic case {(3 = 0) 
to a very strong anharmonicity. Changing the j3 sign the relative phases of 
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Table 6 

Same as table ||, but for the parametric state |<3?i > of eq. (^). In the last column 
the values of the parameter j3 used. 



harm. & lin. 



W 



11 



w 



20 



anharm. 



anharm. & non-lin. 



sin P 



0.26 



0.27 



16.58 



1.96 



1.92 



29.71 



7.20 



-0.02 



0.02 



the |<E>q, > components are changed. The energies of the states were kept fixed 
and equal to the energy, in the harmonic limit, of the main component; i.e. 
Ei = E 2 + + E 3 - and E 2 = E GDR . 

The cross sections corresponding to the |$i > state are shown in table |6| for 
two opposite values of /3. From the table, we can see that the behaviour of the 
cross section is very similar to the one obtained in the complete calculation 
(see table ||). Indeed the results for /3 = —0.02 are similar to the ones obtained 
for the |2 + <g> 3~ > dipole state, where anharmonicities reinforce the effects of 
non-linearities on the cross section. Conversely, for /3 = 0.02 the final result 
is much lower than the one given by the W 20 term alone. This result is very 
similar to the one obtained for the \ISGQR <8> 3~ > state shown in table |5|. 

The reason for this different behaviour can be easily understood in a first order 
calculation if we take into account the following relations 

< u\fi\W(t)\0 > = g E xp(?,t) < u\\V 10 (EX)\0 > 
< M Xfi\W(t)\0 > = , L g EXli (J3,t) < z/ 1 A 1 z/ 2 A 2 |V 20 (EA)|0 >(37) 



1 + 5, 



where g E \p, was defined in equation (p5|). Let us call o\ the cross section 
corresponding to the state |$i >. In a first order calculation we get 



^anharm&non— lin 

<H 

^harm&non— lin 
<*1 



cos P + 



sin/3' 



x 



(3* 



where x is given by the following ratio of matrix elements 

< 2+ ® 3-|V 20 (£l)|0 > 



x 



< GDR\V 10 (E1)\0 > 
Since |x| is usually smaller than 1 the second term in equation 



(39) 
can be 
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important even for small anharmonicities. The values of 8 and x as well as 
their signs are important. 

In the same way we can calculate the cross section er 2 corresponding to the 
state I $ 2 > and get 



^.anharm&non— lin 

— "K 5 t. — = (cos 8 — x sin 8) 2 (40) 

harm&non— lin \ r~ r~ i \ i 

do 



Since \x\ is usually small, the previous ratio will not differ very much from 
one. Note also that, in first order, a ^rm&non-iin and ^harm&iin coinci( ie, while 



a 



harm&non-lin (jjffgj-g from <j narm & :un — Q 



Now, it happens that the x ratio for the |2 + ®3 _ > and \ISGQR®3~ > dipole 
states has the same sign and similar values: —0.058 and —0.091, respectively. 
But the coefficients of their GDR component have opposite sign (see table |2|) 
and their values are such that the dependence of the ratio in equation (|38D on 8 
is nearly linear. Then in one case anharmonicities and non-linearities interfere 
constructively and in the other case interfere destructively. By increasing \8\ 
this property is lost and we could have a reinforce effect in o\ even if (3 and 
x have opposite sign. That is confirmed by the parametric calculation to all 
orders, as can be seen in figure |||, where we have increased the \8\ parameter up 
to about 0.2. That would support that for nuclei with stronger anharmonicities 
we should get a higher increase in the cross section with respect to the linear 
and harmonic case. 

Of course in this example we have assumed that the mixing has just two 
components what is a great simplification, specially for the state whose main 
component is the GDR. Let us go back to the complete calculations where 
we have a mixing of all the states and their proper energies are taken into 
account. 

Similar effects can also be seen on other states in the dipole response of the Pb 
nucleus (see table [|). The two-phonon states located around 25 MeV excitation 
energy are of particular interest. These states are mainly built by coupling the 
giant dipole with the monopole and quadrupole states. As in the previous case, 
the direct transition W 20 and the mixing are important. In addition to that, 
also the transition between one-phonon states contributes to increase the cross 
section. Table ^ shows that in such the increase of the cross section 

of these two-phonon states is more than 300%. This reminds the findings 
discussed in ref ||12|| , where in a very schematic model we were showing that 
non-linearities and anharmonicities might strongly modify the excitation cross- 
section. An example where the anharmonicities play an important role is given 
by the excitation of the \GDR\ ®2 + > state whose cross section is reported in 
table H. The big increase of the anharmonic and non-linear cross section can 
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"o.OO 0.05 0.10 0.15 0.20 
Isin (P)| 



Fig. 6. Relativistic Coulomb target excitation cross section for the parametric 
calculation of eq. ( |36|) as function of absolute value of the mixing coefficient sin (3. 



be entirely ascribed to its big GDR 2 component (see table §). 

This increase of the cross section is seen not only in the dipole channel, which 
gets large contributions from the GDR itself, but also for other multipolarities. 
Let us for instance consider the isoscalar giant quadrupole resonance (ISGQR) 
(see table ||). Looking at its excitation we see that the inclusion of W n raises 
the value of a from 299 mb to 333 mb (in the harmonic case). In this case, 
besides the direct transition to the GQR due to the action of W w , we are 
considering the second order one which proceeds first through the excitation of 
the GDR by W 10 (l) and then to GQR by means of W 11 (see fig. 0). This second 
order process is able to give almost a 12% increase because the excitation 
probability of the first transition is very high and because the effect of W n is 
enhanced by the fact that the energies of GDR and GQR are close each other. 

We close this detailed analysis with a comment on the effect of the non nat- 
ural parity states 1 + and 2~ we have introduced in the calculations. In our 
calculation we found that the contribution of the 1 + states is, in this respect, 
irrelevant and the one of the 2~, in the region of the DGDR, amounts to about 
1 mb. At lower energy, around 16 MeV, its contribution is 2 mb. Similar con- 
clusions have been recently reached in ref . [ 11 ] . 



So far, we have discussed the influence on some particular states. In order to 
get a global view on the effects of both non-linearities and anharmonicities we 
must compute the complete inelastic cross section. Therefore, we have summed 
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GQR 



W "(2) 

g-s. 




b) 



Fig. 7. Schematic representation of the Coulomb excitation of the \GQR > state. 



100.0 - 




Fig. 8. Relativistic Coulomb target excitation cross section for the 208 Pb+ 208 Pb 
system at 641 MeV/A as function of the excitation energy. The three parts corre- 
spond to different energy regions. The cross section for each |3> Q > state has been 
smoothed by a lorentzian with a 3 MeV width. For the low energy region we used 
a 1 MeV width. 
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Table 7 

Comparison between our theoretical results and the experimental cross sections (in 
barn) reported in ref. [6] for the Pb + Pb reaction at 641 MeV per nucleon. The 
theoretical results (first line) correspond to the sum of all GDR (first column) and all 
DGDR (second column) cross-section. The third column contains the cross section 
associated with all the states above the IVGQR (E > 22 MeV). The theoretical 
cross sections are obtained from the non-linear and anharmonic calculation while 
the numbers in parenthesis refer to the linear and harmonic limit. The experimental 
results are reported in the second line. The first number corresponds to the extracted 
GDR cross section while the second number comes from a gaussian fit of the high 
energy cross section after subtraction of the GDR and GQR single-phonon strength. 



GDR 



DGDR 



DGDR energy region 



o"th 



3.13 (3.14) 0.21 (0.22) 



0.31 (0.28) 



<7, 



cxp 



3.28 ±0.05 



0.38 ± 0.04 



up all the contributions coming from the various states after a smoothing of 
each individual line shape by a lorentzian. The results are presented in fig. 
[§]. For the low energy region (fig. |^.a) the width of the lorentzian has been 
chosen equal to 1 MeV, while for the energy region around the GDR (fig. ||b) 
and the one around the DGDR (fig. g.c) it has been fixed equal to 3 MeV. 
In this figure we can see that the single GDR region is not much affected by 
the anharmonicities and non-linearities while the cross-section in the DGDR 
region is increased by 10% when the anharmonicities and non-linearities are 
taken into account. We would like to point out that this increase is mainly 
due to the excitation of two-phonon states whose energies are in the DGDR 
region and whose population has been possible only because of the presence 
of the anharmonicities and the non-linear terms W u and W 20 in the external 
field. The low lying part of the spectrum is also affected and in particular, as 
we discussed before, a new dipole strength is visible in the 9 MeV region. 

In table [7| we show a comparison between our theoretical results and the exper- 
imental cross-section for the GDR and the DGDR energy region. The agree- 
ment for the GDR seems satisfactory. The theoretical yield associated with 
the DGDR states explains about 60% of the experimental cross section. How- 
ever, this disagreement between the experimental cross section in the DGDR 
region and our theoretical estimate is reduced to 18% ± 10% by the inclusion 
of all the different multiphonon states considered in our calculation and lying 
above the IVGQR. 



In conclusion, both the introduction of different two-phonon states and the 
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Table 8 

Same as table [l], but for 40 Ca. 



Phonons 



E(MeV) 



EWSR 



GMRi 


0+ 





18.25 


30 


i\/r z? 


n+ 
U 


u 


OO A T 

zz.4 ( 


04 


GDR\ 


1" 


1 


17.78 


56 


GDR 2 


1" 


1 


22.03 


10 


ISGQR 


2+ 





16.91 


85 


IVGQR 


2+ 


1 


29.53 


26 


3- 


3~ 





4.94 


14 


LEOR 


3- 





9.71 


5 


HEOR 


3 





31.33 
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inclusion of anharmonicities and non-linearities are bringing the theoretical 
prediction rather close to the experimental observation for the Coulomb exci- 
tation of Pb nuclei in the DGDR region. 



5 Results about the excitation of 40 Ca 



We have also done calculations for the excitation of 40 Ca by a 208 Pb projectile 
with Eiab = lOOOMeV/A. The one-phonon basis for 40 Ca are shown in table 
||. We do not have any collective low lying 2 + state because the RPA does 
not generate such state for the 40 Ca nucleus. The properties of the dipole 1~ 
states are reported in table |9], which is the analogous of table |2] for 208 Pb. We 
note that we have bigger anharmonicities than in the Pb case. 

The coupled channel equations (|2~3"1) were solved only for the natural parity 
states which, as we have seen in the case of Pb, are providing the largest 
contribution to the cross-section. The resulting cross section, after a smoothing 
by a lorentzian with a 3 MeV width, is shown in fig. |^. The peak at around 
18 MeV is due to the GDRi with the contribution of the ISGQR state. The 
shoulder at about 22 MeV is given by the GDR 2 and the two-phonon dipole 
state \ISGQR <S> 3~ > which gives, in the anharmonic and non- linear case, a 
10% increase. The latter state is excited in the same fashion of the \ISGQR® 
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Table 9 

Same as table |2|, but for 40 Ca. 



Dipole 



States 



E (MeV) AE (AE 



CGDRx 



cgdr 2 



GDRi 






17.780 


-432. 


0. 


0.989 


-0.006 


GDR 2 






22.034 


-391. 


0. 


0.004 


0.990 


ISGQR 


<8 


3- 


21.851 


708. 


713. 


0.024 


0.011 


ISGQR 




LEOR 


26.616 


231. 


224. 


0.011 


-0.011 


IVGQR 


® 


3~ 


34.541 


-125. 


-128. 


0.001 


-0.020 


GDRi 


<8 


ISGQR 


34.690 


139. 


35. 


-0.063 


-0.044 


GMRi 


<8 


GDRi 


36.026 


-110. 


-214. 


-0.075 


-0.004 


GDR 2 


® 


ISGQR 


38.943 


-21. 


-74. 


-0.034 


0.034 


IVGQR 


<8 


LEOR 


39.305 


-245. 


-245. 


0.000 


0.003 


GMRi 


<8 


GDR 2 


40.280 


-175. 


-292. 


0.011 


-0.079 


GMR 2 


<8 


GDRi 


40.249 


9. 


-202. 


-0.098 


-0.005 


GMR 2 


<8 


GDR 2 


44.502 


20. 


-194. 


0.000 


-0.098 


GDRi 


(8 


IVGQR 


47.379 


-315. 


-308. 


-0.011 


-0.003 


ISGQR 


(8 


HEOR 


48.240 


-13. 


-27. 


0.000 


0.005 


GDR 2 


(8 


IVGQR 


51.633 


-270. 


-271. 


0.001 


0.001 


IVGQR 


(8 


HEOR 


60.929 


-271. 


-275. 


-0.009 


0.005 



3~ > state of Pb, see table [10], with the difference that now the increasing 
factor is 1000 while in the 208 Pb case it was only 100. Finally, we note two 
interesting energy regions where there is a difference between the harmonic 
and linear and the anharmonic and non-linear case, namely the regions around 
35 and 40 MeV. The sum of the cross section for the 1~ state belonging to 
these two regions are reported in table for the different kinds of calculations 
we can make within our approach. From the table we see that the increase is 
essentially due to the dipole 1~ states and this is an almost pure anharmonic 
effect. The global increase in the DGDR energy region amounts to a 20%, 
which is twice what we obtained for Pb. 
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Table 10 

Same as table || but for 40 Ca. 



States 




harm. 




w 20 






& lin. 







anharm. 



anharm. 
& non-lin. 



ISGQR 



l~ 0.004 0.006 6.660 0.284 



3.955 



34 < E < 36 (MeV) 



0.110 



0.287 



0.295 



1.723 



2.221 



38 < E < 45 (MeV) 



0.008 



0.020 



0.022 



1.468 



1.698 




Fig. 9. Relativistic Coulomb target excitation cross section for the Pb+ Ca 
system at 1000 MeV/A as function of the excitation energy. The cross section for 
each |$ a > state has been smoothed by a lorentzian with a 3 MeV width. 

6 Discussion and conclusion 



We have employed an RPA based approach to compute the anharmonicities: 
We have diagonalized the residual interaction between RPA phonons in the 
space of one- and two-phonons states. We have taken into account also the 
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particle-particle and hole-hole terms in the external field making possible the 
direct excitation of two-phonon states as well as the transition between one- 
phonon states. These non-linear terms and the anharmonicities are not taken 
into account in the "standard" approach of the multiphonon picture. Within 
this framework we have calculated the Coulomb excitation of 208 Pb and 40 Ca 
nuclei due to the impinging 208 Pb nucleus at 641 and 1000 MeV/A. 

In this paper we have shown that the inclusion of both anharmonicities and 
non-linear terms in the external field reduce the disagreement between the 
experimental cross section in the DGDR region and the theoretical one calcu- 
lated within the "standard" approach. Moreover, for the 208 Pb case, we have 
found a big effect also at low energy where the |2 + <g) 3 + > state would have 
never been excited without the presence of both the anharmonicities and non- 
linearities. Since these low lying two-phonon states are strongly mixed and 
since their energy is low, we believe that they could also be strongly excited 
by the nuclear part of the mean field at an incident energy lower than the 
one considered here. Theoretical and experimental work in this direction are 
called for. 



In view of our calculations it is clear that non-linearities and anharmonicities 
have an influence on the Coulomb excitation. On some particular states this 
influence can be very strong, while averaging over all the states we have found 
an increase of the cross section by about 10% for Pb and 20% for Ca in the 
region of the two-phonon states while the energies were modified only by a 
few percent. However, this might not be the final answer because of different 
reasons. First, we are working in a truncated subspace in order to keep only 
one- and two- phonon states. However, we know that a large part of the in- 



crease observed in ref. |12| is due to the increase of transition matrix elements 
coming from components of the wave function containing large phonon num- 
bers. These components are not taken into account in the present calculation 
and this reduces the influence of anharmonicities. In fact we have tested this 



point on the simple model reported in reference [I2| and we have observed 
that a truncation of the multiphonon space at the two-phonon level reduces 
the increase of the cross section by almost a factor 4. Unfortunately, this point 
is not easy to improve because the computation time will become too long if 
we are forced to include more multiphonon states. We are now trying to de- 
velop an alternative approach based on time-dependent mean field theory in 
the boson representation. 



The anharmonicities we are computing are mainly due to the residual interac- 
tion into channels which are different from the usual particle-hole interaction. 
One may argue that, as far as effective interactions are concerned, their pa- 
rameters are only fitted close to the ground state. Therefore, except for the 
particle-hole channel the other parts of the interaction are not really con- 
strained by the theory. However, in some cases the residual interaction has 
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been tested far from the ground state. In this respect the relative success of 
the time-dependent mean field theory (and other treatments such as adiabatic 
TDHF or generator coordinate method) might be an indication that the same 
Skyrme parametrisation also holds for large amplitude motion. However, this 
point is certainly calling for more theoretical developments in order to better 
define the effective interaction in channels different from the standard particle- 
hole ones. 

From the experimental point of view it seems that the 208 Pb nucleus behaves as 
a rather good vibrator. In fact the discrepancy about the cross-section between 
theory and experiment is apparently much smaller in the case of Pb than in 
the case of Xe ||. Moreover, as far as the shift in energy of the two-phonon 
state with respect to the harmonic limit is concerned, a shift of less than 
1% was found experimentally for Pb while for Xe it was of the order of 10% 
In our calculation the Pb also appears as a rather good vibrator and the 
predicted effects on the energy shifts are consistent with the experiment. This 
is probably related to the fact that 208 Pb is a double-magic nucleus. It would 
be very important to study non-magic, open-shell or deformed nuclei which 
are expected to be poorer vibrators than double-magic nuclei. In particular, it 
is known that the energy of the GDR is strongly affected by the deformation, 
indicating a possible strong coupling between dipole and quadrupole degrees of 
freedom. This may induce a modification of the cross-section stronger than the 
predicted 10 % to 20% for spherical-magic nuclei. In this respect, extensions 
of the presented results to open shell and deformed nuclei are called for. 

In conclusion, we would like to stress that in addition to the DGDR excitation 
several states are contributing to the cross-section in the DGDR energy re- 
gion. When the non-linearities and the anharmonicities are taken into account 
the total theoretical cross-section above the IVGQR come rather close to the 
experimental result for the Coulomb excitation of Pb. 
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7 Appendix 

We just need to know 
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+oo 



H x ,(/3,t)= / {e^\-l) x ^ + e-^)LO X K^u)dLO (41) 



with \i > 0, since Hxfj.(P,t) = Hx\fi\{P,i). Considering the cases A + /i even 
or odd, together with t positive, negative or null, the integral ( |41| ) will be 
proportional to integrals in |p4]j . Combining all cases we get 



Hx^,t) 



(i + (-i)^)l^r(i±|±H)r(i±^ i )F(i±|±i il i±^ ; i ; -^) - 



i (1 _ (_i)^ ) ^ r(2± ^ )r(2± ^ )F( H±^± ±i ^ ; 3 ; _^) (42) 



These hypergeometric functions can be transformed following [20] as 



F(n + i, n + | — /i; m + ±; — x 2 ) ■ 



(1 _|_ ^2^2n-^-m+l/2 

x F{m — n, m — n + fi; m + ±; — x 2 ) (43) 



If A + fi is even, then m = and n = — Whereas if A + /i is odd, then m = 1 
and n = Therefore, in the cases in which we are interested m — n is an 

integer < 0, and the latter hypergeometric function reduces to the following 
polynomial 
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